Nonlinear modulation of transverse dust lattice waves 
in complex plasma crystals* 
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The occurrence of the modulational instability (MI) in transverse dust lattice waves propagating 
in a one-dimensional dusty plasma crystal is investigated. The amplitude modulation mechanism, 
which is related to the intrinsic nonlinearity of the sheath electric field, is shown to destabilize 
the carrier wave under certain conditions, possibly leading to the formation of localized envelope 
excitations. Explicit expressions for the instability growth rate and threshold are presented and 
discussed. 

PACS numbers: 52.27.Lw, 52.35.Fp, 52.25.Vy 



a. Introduction. Studies of numerous collective pro- 
cesses 01 in dust contaminated plasmas (DP) have been 
of significant interest in connection with linear and non- 
linear waves that are observed in laboratory and space 
plasmas. An issue of particular importance is the forma- 
tion of strongly coupled DP crystals by highly charged 
dust grains, for instance in the sheath region above 
a horizontal negatively biased electrode in experiments 
1^ iSii 4, 5] . Low- frequency oscillations occurring in these 
mesoscopic dust grain quasi-lattices, in both longitudi- 
nal and transverse directions, have been theoretically 
predicted IE 0j !§[ Wi [Tol. a nd later experimentally 
observed [I El iTaTEllllIll ■ We note that the ob- 
servation of the characteristics of transverse vibrations 
around a levitated equilibrium position, where the elec- 
tric and gravity forces are in balance, has been suggested 
as a diagnostic tool, enabling the determination of the 
grain charge 0, 0, 0|- Recent generalizations taking 
into account dust charge variations |18| . layer coupling 
[T^ and two-dimensional crystal anisotropy are also 
worth mentioning. 

It is known from solid state physics that lattice vi- 
brations are inevitably subject to amplitude modulation 
due to intrinsic nonlinearities of the medium. Further- 
more, the wave propagation in crystals are often char- 
acterized by the Benjamin-Feir-type modulational insta- 
bility (MI), a well-known mechanism for the energy lo- 
calization related to the wave propagation in nonlinear 
dispersive media. The MI mechanism has been thor- 
oughly studied in the past, mostly in one-dimensional 
(Id) solid state systems, where nonlinearities of the sub- 
strate potential and/or particle coupling may be seen to 
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destabilize waves and possibly lead to localized excita- 
tions (solitary waves) j22,^3- In the context of plasma 
wave theory, this nonlinear mechanism has been investi- 
gated in a variety of contexts since long ago [23.l25l |. In a 
weakly coupled dusty (or complex) plasma (DPV in par- 
ticular, new electrostatic wave modes arise |ll| 2 6|. vy hose 
modulation has been studied quite recently |27l l28|| : in- 
stability conditions were shown to depend strongly on 
modulation obliqueness, dust concentration and the ion 
temperature pl |. 

In principle, nonlinearity is always present in dusty 
plasma vibrations, due to the form of the inter-grain 
electrostatic interaction potential, which may be of the 
Debye-Hiickel-type H I3 or else [H H^. Further- 
more, the electric potential dominating oscillations in the 
transverse direction, yet often thought to be practically 
parabolic near the levitated equilibrium position jsj, is 
intrinsically anharmonic [s^l , as suggested by experimen- 
tal results 13, "33]. Despite this evidence, knowledge of 
nonlinear mechanisms related to low-frequency DP lat- 
tice modes is still in a preliminary stage. Small am- 
plitude localized longitudinal excitations (described by 
a Korteweg-de Vries equation) were considered in Refs. 

13 based on which longitudinal dust lattice wave 
(LDLW) amplitude modulation was considered in Ref. 
|35| . However, to the best of our knowledge, no study has 
been carried out, from first principles, of the amplitude 
modulation of transverse dust lattice waves (TDLWs) be- 
cause of the sheath electric field nonlinearity. This Letter 
aims in making a first analytical step towards the study 
of DP crystals in this framework. 

About three years ago, Misawa et al. ^3 reported the 
observations of vertical nonlinear oscillations of a dust 
grain in a plasma sheath, and interpreted their results 
in terms of a position-dependent delayed grain charg- 
ing effect. Zafiu et al. ^3!| studied vertical dust-grain 
oscillations in a rf-discharge plasma, and succeeded in 
pointing out their strongly nonlinear behavior due to the 
sheath potential anharmonicity. Recently, Ivlev et al. 
36] investigated the nonlinear coupling between high- 
frequency transverse (vertical) dust lattice oscillations 
(TDLOs) and slow longitudinal dust lattice vibrations 
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(LDLVs). Within the framework of a slowly varying en- 
velope approximation, they derived a pair of equations 
for the modulated TDLOs and driven (by the pondero- 
motive force of the latter) slow LDLVs. The coupled 
system of equations admit envelope soliton solutions. Fi- 
nally, compressional pulses were studied, both theoret- 
ically and experimentally, in Ref. 37]; however, these 
excitations are related to longitudinal grain motion. 

In this Brief Communication, we consider the am- 
plitude modulation of transverse dust lattice waves 
(TDLWs), taking into account self- interaction nonlin- 
earities associated with harmonic generation, involving 
the intrinsic nonlinearity of the sheath electric poten- 
tial. Thus, the underlying physics of our nonlinear pro- 
cess is entirely different from those considered in Refs. 
[13. l33. j37j . By adopting the reductive perturbation 
technique, we derive a cubic nonlinear Schrodinger equa- 
tion for the modulated TDLWs. It is shown that the 
latter may be modulationally unstable depending on the 
plasma parameters, and that they can propagate in the 
form of envelope localized excitations due to a balance 
between nonlinearity and dispersion. Explicit forms of 
localized excitations are presented. 

b. Equation of motion. Let us consider TDLWs 
(vertical, off-plane) propagating in the one-dimensional 
(Id) DP crystal. The dust grain charge q and the mass 
M are assumed to be constant. DP crystals have been 
shown to support low-frequency optical-mode-like oscil- 
lations in both transverse and longitudinal directions 
Q IE IE IE IE ^1 ■ Focusing on the former and summariz- 
ing previous results, let us recall that transverse motion 
of a charged dust grain (mass M, charge q, both assumed 
constant for simplicity) in a DP crystal (lattice constant 
rg) obeys an equation of the form 



M- 



= ^'^ ^0 (2 - 6Zn-l - 6Zn+l) + F, 



Mg, 
(1) 

where 6zn = z„ — zq denotes the small displacement of 
the n— th grain around the equilibrium position zq, in the 
transverse (z—) direction, propagating in the longitudi- 
nal (x—) direction. Assuming that the neighboring dust 
grains (situated at a distance \) interact via 

an electrostatic potential ^(x), we obtain the DP oscilla- 



tion 'eigenfrequency' cJq = — 



Mro dx 



e.g. m 



the 



case of a Debye-Hiickel potential: $(x) — [q/x) e 



— x/\d 



UJn 



12. ] p-ro/Xn 



(2) 



where A d denotes the effective DP Debye radius 0| . The 
force Fe = qE{z) is due to the electric field E{z) = 
—dV{z)/dz; the potential V{z) is obtained by solving 
Poisson's equation, taking into account the sheath po- 
tential and also (in a more sophisticated description) the 
wake potential generated by supersonic ion flows towards 
the electrode |T3|. The potential V{z) thus obtained, ac- 
tually a nonlinear function of z, can be developed around 



the equilibrium position zq as 

V{z) « V{zo) + V(i) Sz + ^ ^(2) {6zf + i ^(3) {Szf 

+0[iSzr]; (3) 



the electric force therefore 



obviously V^j) 
reads 

Fe(z) « Fe(zo)+7(i) 5z+^(2) (<5z)'+7(3) {5zf+0[{5zf] , 

where all coefficients are defined via the derivatives of 
V{z), i.e. 7(j) = —q\V^j^iy The zeroth-order term bal- 
ances gravity at (and actually defines the value of) Zq, viz. 
Feizo)-Mg = 0, while -7(1) = qV^2) = 1 = M uj^ is the 
effective width of the potential well; the value of the gap 
frequency cug = uj{k = 0) may either be evaluated from 
first theoretical principles |E| or determined experimen- 
tally ^E|j and is typically of the order of ujg/2TT w 20 Hz 
in laboratory experiments. Collisions with neutrals and 
dust charge dynamics are omitted, at a first step, in this 
simplified model. 

Retaining only the linear contribution and considering 
phonons of the type, Xn = An exp[i (knrQ — ujt)] + c.c, 
we obtain an optical-mode-like dispersion relation 
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(4) 



where uj and k — 2t:/X denote, respectively, the wave fre- 
quency and the wavenumber. We will not go into further 
details concerning' the linear regime, since it is sufficiently 
covered in Refs. (E [tI, |8|, |9|, |lE| ■ Let us now see what hap- 
pens if the nonlinear terms are retained. 

c. Derivation of a Nonlinear Schrodinger Equation. 
For analytical tractability, we shall limit ourselves to 
a quasi-continuum limit, by considering an amplitude 
which varies over a scale L which is significantly larger 
than the inter-grain distance tq (i.e. L/rQ <^ 1). Equa- 
tion takes the form 



d^Ur, 



-I-Cq {2 5Un- Un-l 



U„+l)+CJg U„ + Q!?4+/3ufj = 0, 

(5) 

where we set 5z = u{x,t) for simplicity; cq = ujq tq is 
a characteristic propagation speed related to the interac- 
tion [e.g. Debye-type, see Q] potential; the nonlinearity 
coefficients a, (3 are related to the anharmonicity of the 
electric potential, viz. 

« = -^-^, /3 = -?#-^- (6) 



M 2M 



M 6M 



Remember that inter-grain interactions are repulsive, 
hence the difference in structure from the nonlinear 
Klein-Gordon equation used to describe one-dimensional 
oscillator chains. 'Phonons' in this chain are stable only 
in the presence of the electric field (i.e. for 7 7^ 0). 
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We now proceed by considering small-amplitude oscil- 
lations of the form 



u — eui + u-i' -I- ... 



at each lattice site. Introducing multiple scales in time 
and space, i.e. = e"a;, T„ = e"t (n = 0, 1, 2, 

...), we develop the derivatives in Eq. (jSJ in powers of 
the smallness parameter e and then collect terms aris- 
ing in successive orders. The equation thus obtained in 
each order can be solved and substituted to the subse- 
quent order, and so forth. This reductive perturbation 
technique is a standard procedure for the study of the 
nonlinear wave propagation (e.g. in hydrodynamics, in 
nonlinear optics, etc.) often used in the description of lo- 
calized pulse propagation, prediction of instabilities, etc. 
[2^ I23I I . This procedure leads to a solution of the form 



t) 



2 1^1^ 



+ (7) 



where c.c. denotes the complex conjugate; recall that lo 
obeys the dispersion relation Q). 

The slowly- varying amplitude A — A{Xi — VgTi) 
moves at the (negative) group velocity Vg — diu/dk = 
— Wq^o sin(fcro/ci;) in the direction opposite to the phase 
velocity; this backward wave has been observed experi- 
mentally: see the discussion in Ref. 14]. The amplitude 
A obeys a Nonlinear Schrodinger Equation (NLSE) of 
the form 



dA 
df 



P 



d^A 



OlApA^O, 



(8) 



where the 'slow' variables {X,T} are {Xi — VgTi,T2}, 
respectively. The dispersion coefficient P, which is re- 
lated to the curvature of the phonon dispersion relation 
Q as P = (d2w/dfc2)/2, reads 



P : 



9 9 

4w3 



2 I cos(fc7-o) +cos(2A:ro) 



and the nonlinearity coefficient 



Q 



1 /lOa^ 



2u; V 3^2 



3/3 



1 



2Muj 



(2) 



37(1 



3 7(3) 



(9) 



(10) 



is related to the electric field nonlinearity considered 
above. Notice that the sign of P depends on the ra- 
tio Wg/wp = A > 0, which appears naturally as an 
order parameter; see for instance in Q that A > 2 
is a stability criterion (necessary for lo to be real in 
the whole range of the first Brillouin zone). For long 
wavelength values, P « —{clLjg)/{2ui^) < 0, given the 
parabolic form of uj{k) close to fc = (continuum case). 
In the general (discrete) case, we see that the coeffi- 
cient P becomes positive at some critical value of k, say 



kcr, inside the first Brillouin zone. Some simple algebra 
shows that the zero-dispcrsion point kcr satisfies the re- 
lation: cos{kcrro) = (2 - A^ + A^A^ - 4)/2, for A > 2 
(0 < kcrfQ < tt); otherwise, for A < 2, P remains nega- 
tive everywhere. 

d. Modulational instability. In a generic manner, a 
modulated wave whose amplitude obeys the NLS equa- 
tion ©, is stable (unstable) to perturbations if the prod- 
uct PQ is negative (positive). To see this, one may first 
check that the NLSE accepts the monochromatic solution 
(Stokes' wave) A{X,T) = Aoe*<3l^°l'^ + c.c. The stan- 
dard (linear) stability analysis then shows that a linear 
modulation with the frequency SI and the wavenumber k 
obeys the dispersion relation 



(11) 



which exhibits a purely growing mode for k > K^r = 
{Q / PY^'^ \Aq\. The growth rate attains a maximum 
value of ^max — Q\Aq\'^ . This mechanism is known as 
the Benjamin-Feir instability For PQ < 0, the wave 
is modulationally stable, as evident from Hll|l . 

One now needs to deduce the sign of Q, given by H10(l . 
in order to determine the stability profile of the TDL 
oscillations. In fact, given the above definitions of the 
parameters u, a, one easily finds that Q is related to 
the (derivatives of the) electric potential V{z) via 



= 



q 



5^(3) 



3Vr 



V, 



(2) 



(4) 



3^(2) 



(4) 



(2) 



(12) 



The exact form of the potential V{z) may be obtained 
from ah initio calculations or by experimental data fit- 
ting. For instance, in Ref. the dust grain potential 
energy U{z) = q V{z) was reconstructed from experimen- 
tal data as 



0.9 6z + i {Szf - io.5((5z)3 



+ io.07(fc)4 + 
4 



, (13) 



which is Eq. (9) in Ref. [H; upon simple inspection 
from ©, we obtain V(3)/V(2) = -1 , V(4)/V(2) = 0.42, 
so the value of Q is positive, as may be checked from 
(|12|l . Therefore, the transverse oscillation considered in 
Ref. would propagate as a stable wave, for large 

wavelength values A. However, for shorter wavelengths, 
the coefficient P = uj"{k)/2 - as defined in - may 
become positive (and so will the product PQ, in this 
case); the TDL wave may thus be potentially unstable. 
These results may a priori be checked experimentally. 

e. Localized excitations. A final comment concerns 
the possibility of the existence of localized excitations 
related to transverse dust-lattice waves. It is known 
that the NLSE © supports pulse-shaped localized so- 
lutions (envelope solitons) of the bright {PQ > 0) or 
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dark/grey (PQ < 0) type |3l,|3^. The former {contin- 
uum breathers) are 

A = (2L'/PQ)i/2 sech[{2D/PQy/^ {X - T)] 

X exp[iVe{X -VcT)/2P] + c.c. , (14) 

where Ve (vc) is the envelope (carrier) velocity and D — 
{v1 — 2ve Uc)/(4P^); they may occur and propagate in the 
lattice if a sufficiently short wavelength is chosen, so that 
the product PQ is positive. We note that the pulse width 
L and the amplitude p satisfy Lp ^ (|P/Q|)^/^ = const. 
For PQ < 0, we have the grey envelope soliton [s^ 

A = piil-a^ sech'^{[X-{v^+2a)T]/Li}Y/'^ exp(icr), 

(15) 

where tr — a{X,T) is a nonlinear phase correction to be 
determined. This excitation represents a localized region 
of negative wave density (a void), with finite amplitude 
(1 — a)pi at X = 0; < a < 1). Again, the pulse 
width Li = {\P/Q\)^/^ /{api) is inversely proportional to 
the amplitude pi. Notice the (dimensionless) parameter 
a, which regulates the depth of the excitation. For a = 1, 
one obtains a dark envelope soliton, which describes a lo- 
calized density hole, characterized by a vanishing ampli- 
tude at C = 0. The latter excitations (of grey/dark type), 
yet apparently privileged in the continuum limit (where 
PQ < 0), are rather physically irrelevant in our (infinite 
chain) model, since they correspond to an infinite energy 
stored in the lattice. Nevertheless, their existence locally 
in a finite-sized chain may be considered (and possibly 
confirmed) either numerically or experimentally. 



In conclusion, we have shown that the modulational 
instability is, in principle, possible for transverse DP lat- 
tice waves. Long wavelength modes seem to ensure wave 
stability, while shorter wavelength modes may be mod- 
ulationally unstable. The existence of localized excita- 
tions and the occurrence of the modulational instabil- 
ity rely on the same criterion, which needs to be thor- 
oughly examined for a given exact form of the sheath 
electric potential. These results may be investigated and 
will hopefully be confirmed by appropriate experiments. 
Of course, for a more complete description, one has to 
take into account certain factors that are ignored in this 
simple model, viz. collisions with neutral particles, dust 
charge variations and, eventually, transverse to longitu- 
dinal mode coupling. The effect of inter-layer coupling 
may also be investigated. Work in this direction is in 
progress and the results will be reported soon. 
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